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ABSTRACT 

We present a "feedback compression" model to describe the galactic spheroid 
formation and its relation with the central nuclear activity. We suggest that the 
star formation itself can serve as the "positive feedback" in some extremely dense 
region to trigger the starburst. The star formation rate as well as the related 
stellar feedback-induced turbulence will be maximized under the regulation of 
the background dark halo's gravity. There is also stellar feedback acting inward 
to confine and obscure the central black hole (BH) till the BH grows sufficiently 
large to satisfy a balance condition between the accretion disk wind and the 
inward stellar feedback. The extremely vigorous star formation activity, the BH- 
bulge relations, the maximum velocity dispersion as well as the maximum BH 
mass are investigated based on such scenario, and are found to be consistent with 
observations. 

Subject headings: black hole physics - galaxies: formation - galaxies: nuclei - 
galaxies: starburst - galaxies: structure 



INTRODUCTION 



Observations have shown clear evidence that the mass of supermassive black hole (SMBH) 
in the center of every galaxy is tightly correlated with the velocity dispersio n of the bulge 



stars ( Ferrarese fc Merritt 2000:lGebhardt et al. 



2000 



Tremaine et al. 2002) and the mass 



of the whole bulge (iMagorrian et al 



1998 



iMcLure fc Dunlop 



2002 



Marconi fc HuntI 120031 ) 



which reflect the interplay with the nuclear activity and the spheroidal star formation ac- 
tivity (potential well depth of bulge stars). Recent consensus emphasizes that the nuclear 
black hole (BH) feedback is the sticking point to explain these correlations: the central BH 
can interact with the surrounding environment through BH feedback in a self-regulated way. 



- 2 - 



However, the tightness of the correlations indicated by both the observations and numerical 
simulations pose big challenge to many conventional feedback models, whose results are sen- 
sitive to the variable par ameters such as gas fraction (the mass fraction of gas to dark matter) 
( iDi Matteo et al. II2005I ). Recently, some mod i fied fe edback models have been proposed to 
alleviate the inconsistency. iBegelman fc Nath I (120051 ) suggest that the BH accretion physics 
on local scale (relative to previous "global" models) should be considered. They consider 
that the " maximized" accret ed gas mass will make the result insensitive to the gas frac- 



tion. IXu. Wu fc Zhao I (120071 . hereafter XWZ07) provide another scenario to explain both 
the starburst activity in protospheroid and the BH-bulge relations. They suggest that the 
spheroid formation is "halo-regulated": the star formation is regulated by the background 
dark matter halo. In particular, the strong star-forming feedback acts outward to resist 
the gravity from the dark halo while acts inward to feed and obscure the central BH. Such 
scenario can naturally link the nuclear activity and outside star formation activity, and well 
explain the BH-bulge relations. 

Here we revisit the idea of XWZ07 and find the " halo- regulated" mechanism can be 
generalized and simplified. We argue that it is possible that the stellar feedback can compress 
the surrounding gas and trigger the catastrophic star formation activity in a very dense 
environment when neglecting the heating effects. The intense starburst will make the star- 
forming region highly turbulent. The "maximized" turbulent sound speed, which is regarded 
as a measure of potential well depth of the protospheroid, can be reached when the whole 
system is in the virial equilibrium. Based on such "feedback compression" model, we show 
that once the assumption of homogeneous turbulent environment is valid, the star formation 
is regulated by the dark halo's gravity and regulates the central BH growth. The resultant 
BH-bulge relations are universal: only related to the velocity dispersion of stars and dark 
matter profile. 



2. STAR FORMATION UNDER THE MOMENTUM FEEDBACK 

COMPRESSION 



The merger induced starburst regions h ave very dense environmen t: the molecular gas 
density can be as high as 10^ ~ lO'^cm"^ (IDownes fc Solomon Ill998l ). much denser than 
that of disk galaxies. The cooling by collisionally excited atomic and molecular emission 
processes can be very efficient. Hence, we hypothesize that the heating feedback due to 
the ionization of HII region may be neglected comparing with the mechanical feedback in 
the starburst region. Since the shock-heated gas can quickly radiate their thermal energy, 
it is possible that the outward propagating shock can compress the surrounding gas into a 
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gas shell in a momentum driven way and trigger the new star formation events. The result 
can be understood by a simple example: assuming that the expanding gas shell is thin, its 
dynamics can be described by 



dt V dt J Airpo ' ^ ' 

where is the mass loss rate, is the velocity of wind and po is the density of homogeneous 
gas around the stars. We can easily find the following solution at large R 




Kit) = t'^'. (2) 



The shell will finally stall due to the ambient pressure, and we have the stalling radius 
Rs = {Mn,v^/AnnokTY^^. Using Eq. (2), the total compression timescale can be estimated 

as 

/ • \ 



4fcT \ Snp. 







Comparing with the dynamical timescale t^yn = (37r/32Gp)^'^, we have 

-sf^y'A'e^fc', (4) 



tdyn V Po 

where pc is the gas density of the compressed gas shell, Mq = Mw/lO^^MQ-yr^^, v^qq = 
Vw/500kms~^ and Tiq = T/IOK. Note that the shell gas density pc is larger than pq, so tc 
is at least one order larger than tdyn- It demonstrates that in a very dense environment, the 
density inhomogeneities are amplified by the feedback compression and the gas has enough 
time to collapse into the cloud. 

The " light" outflow accelerating against the dense cold gas shell will trigger the Rayleigh- 
Taylor instability (RTI). Moreover, the Kelvin-Helmholtz instability (KHI) may also occur 
as the outflow punches into the gas shell and moves through the cold gas. The combination 
of these two instabilities will prevent the gas collapse and destruct the star-forming cloud. 
If the densit y contrast is large, the dispersion relation for RTI is w; ~ (/ca)^/^ where a is the 



acceleration (iChandrasekhar II1961I ). The characteristic growth timescale for RTI responsible 
for the cloud destruction is trt ~ (2vr/ j/a)^/^ where Ij is the Jeans length. The acceleration 
can be derived from Eq. (2) 

4 \ 2 vrpo / 
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Initially, the growth timescale trt may be smaller than the dynamical timescale tdyn- The 
RTI will make the dense cold shell highly porous to the hot w ind and entrain cold gas into 
the wind. The efficiency of driving wind enhances (ISilk Il200ll ). Note we define 

the duration in which RTI dominates by setting trt = idyn- Using Eq. (5), the duration is 

^1/3 

yr, (6) 



2.5 X lQ^Ml'\ll,W%, 



dyn- 
1/3 I Pc 



Po 



where = hq / 10^ cm~~^. The duration is smaller than the dynamical timescale tdyn, which 
indicates that t he cloud destr uction by RTI and KHI (its growth timescale is of the same 
order as RTFs (lAgertz 1 12006| )) may not be important in the momentum driven phase, at 
least at large radius R. So it is reasonable to assume that the feedback induced compression 
can trigger and expedite the cloud collapse and the star formation. 



"HALO-REGULATED" STAR FORMATION AND BLACK HOLE 

GROWTH 



Star formation in the local disk galaxies is us ually inefficient and the typical star for- 
mation efficiency (SFE) is ~ 2% (IKennicutt Ill998l ). It is because that the star formation is 
regulated by the "negative" feedback (heating, blown out) to keep the star formation rate 
(SFR) from raising too high. However, as mentioned in the above section, it may not be the 
case in the starburst region. Without the heating effects, the feedback induced compression 
can serve as a "positive" feedback to trigger the intense star formation activities unless the 
total mechanical energy generated by the stellar feedback is larger than the binding energy 
of these regions. In another words, once the total feedback from those massive stars are 
unable to disrupt or unbind the whole star forming region, the star formation process may 
continue for a relatively long time, which is analogous to the formation of bound clusters 
( lElmegreen &: Efremov 1119971 ). 



Following XWZ07, we adopt the NFW density profile for the background dark matter 
and assume the standard cosmological parameters with Qm = 0.3, Q\ = 0.7 and h = 0.7 
{z = 0). The inner NFW profile is given by 



pNFw[r) 



n = ISOMgpc 



c,0.58 



(7) 



flNavarro. Frenk fc White 1119971 : iKomatsu fc Seliakll200l[ ). where M^^u 



is the virial mass of the halo, ^{z) = \(flm/^mi z))(ArJlOO)], fl 



Ac = 187r2+82d-39rf2 where d = ^^(2)-! flBrvan fc Norman lll998l : lBarkana fc Loeb 



M^ir/IO^^MQ and 



200lh . ^ 



c,0.58 



[ln(l + c) - c/(l + c)]-V 0.58 where c 



13.4M-°2"( 



1 + z) ^ is the con- 
centration parameter (IBuUock et al. II2OOII ). Such profile in the inner region is almost a 
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universal profile: it doesn't change by the galaxy merger or interaction, it is also insensitive 
to the redshift and the halo mass. Another interesting thing to note is that such inner profile 
gives a nearly constant gravitational acceleration 

g^^(r) = ^^l^f ^""^ = 2nGU ~ 1.2 x IQ-^cm sec'^. (8) 

Intense star formation activity and related stellar feedback will make the protospheroid 
environment highly turbulent. Because the bulg e stars density distribu tion implies that the 



protospheroid density profile may follow p ~ r ( iTremaine et al. Ill994l ). here we assume an 



isothermal density profile for the protospheroid baryon or mixture of gas and stars 

where pb{r) is the total baryon density or density of the mixture of gas and stars, Cg is the 
"isothermal" turbulent sound speed and Mj, is the enclosed mass within the radius r. The 
momentum transport during the compression requires that 

Pb{r)cl, (10) 



47rr^ 



where is the net outward momentum deposition rate. The star formation activity and 
the feedback-induced turbulence will be maximized till the whole system evolves to a virial 
equilibrium state. Once the whole system becomes a little more turbulent, the deposited 
momentum flux will make the whole system deviate from the equilibrium state and the SFR 
is hence avoided from raising higher. Such a self-regulated mechanism will make the whole 
system maintain the equilibrium state. 

We can write the equation of the virial equilibrium for the protospheroid as 

GAf'^ 

SMbC^ ^ = + irUGMbT, (11) 

r 

where Cs^m is the maximum turbulent sound speed. The first term of the right side of Eq. 
(11) denotes the total baryon's self-gravity while the second term denotes the gravity from 
the background dark matter. Substituting Eq. (9) into Eqs. (10) and (11), we obtain 

A = = nllGMh. (12) 

Eq. (12) shows that the maximum turbulent sound speed or the potential well depth of the 
protospheroid is directly related to the background dark matter. In another words, both the 
star formation activity and the feedback compression are regulated by the dark halo's gravity 
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as XWZ07 proposed. It only requires the homogeneous turbulent environment and the virial 
equilibrium, without involving the detailed gas assembly physics such as monolithic collapse 
or merger driven inflow. So the feedback compression combined with the halo regulation 
scenario provides a more general description to the star formation process during the spheroid 
formation. 

Following XWZ07, during the formation of the protospheroid, the stellar feedback can 
act in both inward and outward directions. At small scale (e.g. galactic nuclear region), 
the inward stellar feedback (required to conserve the local momentum) obscures and regu- 
lates the BH growth, while the outward stellar feedback resists the gravity at large scale. 
In particular, the inward stellar feedback regulates the BH growth by interacting with the 
Compt on-thick wind launched from the accretion disk if super-Eddington accretion is as- 



sumed ( iKing &: Pounds II2003I ). The final balance between the inward stellar feedback and 



the disk wind is achieved when 

p. = (13) 

K 

where is the momentum flux transported by inward stellar feedback and the right side of 
Eq. (13) is the momentum deposition of the disk wind. Then if the BH's feedback is large 
enough to halt the further gas supply, and a BH will end its main growth phase after Eq. 
(13) is satisfied. 

At a late epoch of the coeval evolution, star formation consumes most of gas and 
gradually fades away. Without continuous ejecting energy and momentum from the stel- 
l ar feedback, the fe edback-induced turbulence will decay on a crossing time of the system 



(jStone et al. Ill998l ). Then the remaining stellar system will be virialized through violent 
relaxation under the combined gravity from itself and the background dark matter. We call 
the system at "initial state" after turbulence decay and before virialization, and at "final 
state" after virialization. Assuming the total energy of the initial state is E, the kinetic 
energy of the final virialized state is i^' = —E according to the virial theorem. So we have 



Using Eq. (12), we obtain 



= 2 ( ^ + vrHGr ) = 6c^,^. (14) 



4 

P. = ^, r,= (15) 

where rj, is the boundary radius of the initial state. We note that the total momentum 
deposition rate from stars is only related to the velocity dispersion of the final stellar system, 
independent of any parameters of the detailed star formation physics. 



- 7- 



Using Eqs. (13) and (15), we obtain the final BH mass 



final 
BH 



167rG2 



1.5 X IOVIooMq. 



(16) 



The result is remarkably consistent with the low-redshift observations (ITremaine et al. 



20021 ) ■ The stellar bulge mass is approximately equal to M^. Using Eqs. (12) and (13), 



the ratio of BH mass to bulge mass can be expressed as 

Mbh f^n 



M, 



bulge 



1.4 X 10-3M°fJ[e(^)]^/^vl/-i 



58! 



(17) 



which matches the Magorrian relation found fo r the local galaxies (iMagorrian et al.lll998 



McLure fc Dunlop II2002I : Marconi fc HuntI 120031 ) . 



4. APPLICATION TO THE HIGH-REDSHIFT STAR FORMING 

GALAXIES 

The gas fraction of high-redshift galaxies is much larger than that of local galaxies. In 
an extremely case, we take the fraction ~ 1. In another words, the protospheroid with mass 
of Mb is almost totally in the gaseous form. Large amount of gas accumulating in the central 
region will trigger the central vigorous starburst accompanying with strong star-forming 
feedback. We mainly focus on two primary sources of star forming feedback: radiation 
pressure and supernovae. The combined momentum flux deposited in these star-formin g 
feedback can be written as (IMurray. Quataert fc Thompson Il2005t IXu. Wu fc Zhao 1120071 ) 



P = P + P 



(1^ 



where M^ is the star formation rate and = 1 + Pgn/Prp is of the order of unity in our 
model. Combining Eqs. (13), (16) and (18), we can easily obtain the star formation rate of 
high redshift starburst galaxies 



(19) 



Although the star formation law in the protospheroid is far from clear, we ado pt an equivalent 



Schmidt-Kennicu tt Law in order to compare with the local disk galaxies ( jSchmidt I Il959 
Kennicutt 1^9981 ) 

Mb , . 

riT^, (20) 



yn 
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where rj is the equivalent star formation efficiency and tdyn — (37r/32G'pg)^/^ is tlie dynamical 
time scale. In our model we take Pg = 3Mft/47rrf as the average gas density where is the 
outer boundary radius given in Eq. (15). 

Prom Eqs. (12) and (18), we have 

CmeM^c = TTUGMh. (21) 

Using Eq. (20) to eliminate M^/M^ in Eq. (21), we get the equivalent star formation 
efficiency as 

V = ^-OAa,oo^-h^\ (22) 

where 63 = e/10~^. 

We find that our derived equivalent SEE is much higher than that in normal disk galaxies 
inferred from the Kennicutt Law, but is consistent with the high redshift star formation 
observations and some small scale star formation observation (eg. SEE in the formation 
of the protocluster) . We note that the derived SEE is independent on the radius of the 
star-forming region and time, and the larger the velocity dispersion is, the higher the star 
formation efficiency is. Eq. (22) also gives us another implication for the maximum stellar 
velocity dispersion. We can rewrite Eq. (22) as 

af = . (23) 

TT 

The physical limit requires 77 < 1, so the maximum stellar velocity dispersion is 



4^26 



^max ^ ^ MO^^€skmS-\ (24) 



TT 



According to Eq. (16), the maximum BH mass is 8 x 10 Mq. 

Using Eqs. (12), (15) and the expression of tdyn , we can also obtain the characteristic 
star formation timescale 



tdyn 2r6^nT,eC ^m^C 

~q' ~ a) ~ ttTIG 
^ 1.0 X IQ^UezM^^ifUzT'^'^-^co.myr. (25) 

Through the feedback compression, the dark halo's gravity regulates the SER and SEE 
to a higher level during the spheroid formation. We note that for some luminous elliptical 
galaxies whose velocity dispersion cr ~ 300A;ms~^, the predicted star formation rate can reach 
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as high a s 3000M(Ttyr~^ , which is consiste nt with the observations of high-redshift starburst 
galaxies (ISolomon fc Vanden Bout II2005I ). It is also interesting to note that the star forma- 
tion time scale is independent with velocity dispersion a. Furthermore, the characteristic 
star formation time scale is larger than the salpeter time scale, which is usually taken to be 
the typical time scale for BH's growth. This means that the BH will first grow relatively 
fast to reach the balance condition (Eq. (13)) and then grow relatively slow to response the 
outside stellar feedback. 



5. DISCUSSIONS 



Unlike the previous momen tum feedback models ( Kingj[2003 : Murray. Quataert fc Thompson 
2OO5I : iBegelman fc Nathlboosh . which mainly focus on the BH feedback dynamics, our sce- 
nario considers more about the star formation activity in the protospheroid and its relation 
with the nuclear BH growth. Our model favors gas-rich environment at high-redshift because 
large amount of gas is needed to form stars and obscure the BH. We argue that at early 
epoch of BH growth (when BH is relatively small), large scale outflow or jet may not be 
crucial in pro ducing BH-bulge relati ons although they become important after BH's main 
growth epoch (jChurazov et al. II2002I ). In another words, BH doesn't generate outflow or jet 
till the balance condition Eq. (13) is satisfied. The derived Mbh — relation based on our 
model is insensitive to the gas fraction and other variables, because the velocity dispersion is 
the measure of the maximized turbulent velocity of the total baryon component rather than 
the gas only, as Eq. (9) shows. The derived Mbh — M^uige relation has weak dependences 
on the redshift and the halo mass, which offer the intrinsic scatters to the relation. 

Extremely high SFR and SEE are the results of certain "posit ive" feedback which is 
probably due to either "internal" or "external" effects. ISilk I (120051 ) suggest that the high 
SFR and SFE are triggered by the super-Eddington outflow driven by the SMBH. Such 
"external" positive feedback naturally leads to a top-heavy initial mass function (IMF) 
which is preferred by the e arly generation of star formation and predicts an antihierachi- 
cal trend of SMBH growth (IMerloni. Rudnick fc Di Matteo 1120041 ). However, recent optical, 
infrared and X-ray studies o f SMGs indicate that SMGs harbour relatively smaller SMBH 



than that of typical qu asars (jivison et al.lll998l : IVernet fc Cimatti II2OOII : ISmail et al. 112003 



Alexander et al. II2005I ) and the main growth phase of SMBH ("pre-quasar" phase) is heavily 
obscured. Considering the SM Gs themselves are massive galaxies which are reckoned as the 
progenitors of local ellipticals (iGreve et al.ll2005l ). the vigorous star formation activity can 
not be regulated by the small BH. Reversely, the star formation activity should have great 
impact on the small BH. So we argue that the "internal" positive feedback which is produced 
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by the star formation itself may exist in some extremely dense starburst region at median 
redshift, and the SMBH outflow t riggered star formation mode may only be available at very 
high redshift (IWalter et al. II2004I ). In addition, the maximized "positive" stellar feedback is 
actually related to the dark halo's gravity in our model. Under the regulation of dark halo's 
gravity, the maximized velocity dispersion also has its maximum value, which is determined 
by the physical upper-limit of SFE. The physical upper-limit of BH mass (~ lO^^M©) can 
then be obtained by the Mbh — o relation , and the dynamical signature of such SMBHs 

2OO3I1 nee do 



should be detectable (IWvithe fc Loeb 



support such result (INetzer I l2003l : IVestergaard II2004I ). although all of them still contain a 
lot of uncertainties. We expect more accurate SMBH mass measurements in the future to 
confirm our result. 
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